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Abstract. - The influence of Gaussian laser pulses on the transport through molecular 
wires is investigated within a tight-binding model for spinless electrons including correlation. 
Motivated by the phenomenon of coherent destruction of tunneling for monochromatic laser 
fields, situations are studied in which the maximum amplitude of the electric field fulfills the 
conditions for the destructive quantum effect. It is shown that, as for monochromatic laser 
pulses, the average current through the wire can be suppressed. For parameters of the model, 
which do not show a net current without any optical field, a Gaussian laser pulse can establish a 
temporary current. In addition, the effect of electron correlation on the current is investigated. 



Introduction. - In recent years many groups have been working on making the vision 
of molecular electronics reality [1,2]. This bottom-up approach for electronic devices has 
certain advantages over standard top-down approaches which are mainly used these days. In 
molecular electronics the transport is through single molecules or molecular aggregates and has 
therefore to be treated quantum mechanically. This quantum nature makes it of course more 
complicated to determine, for example, the current-voltage characteristics than in classical 
theories. But at the same time a quantum treatment offers certain advantages and especially 
the possibility of constructive or destructive interference effects. 

In the current letter we focus on the electron transport through a molecular wire which is 
coupled to two leads acting as electron source and drain. Many theories of quantum transport 
utilize the non-equilibrium Green's function approach [3] which is formally exact within the 
lead-wire coupling. But because of the dependence of Green's functions on two time argu- 
ments it is rather difficult to determine the current without any further approximations, as 
for example, the wide-band limit. Another possible route is to treat the lead-wire coupling 
perturbatively and derive quantum master equations for the electron dynamics in the wire and 
the current through the wire [4-7]. In addition to the wire-lead coupling a laser field can be 
coupled to the wire and/or the leads. This would possibly allow for an ultrafast opto-electronic 
device. First experimental [8] and theoretical [4, 7] investigations in this direction have been 
performed. Most of the theoretical studies are based on a master equation approach since in 
this formalism allows easily to include time-dependent laser fields. For time-periodic fields the 
Floquet theory can be employed [4]. Recently the present authors derived a formalism which 
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Fig. 1 - Sketch of the two-site system defined in the text. 

treats the laser- wire interaction exactly for an arbitrary time-dependence of the field [7] . This 
approach has been developed in analogy to the treatments in dissipative quantum dynamics, 
but instead of a bosonic heat bath in dissipation theories the leads are treated as fermionic 
reservoirs. In addition, the system-bath coupling operator allows for particle exchange in the 
case of molecular wires while in dissipative theories only energy exchange is allowed. The 
present approach makes use of the spectral decomposition of the spectral density J(uS) first 
introduced by Meier and Tannor [9] . These authors applied a time- nonlocal formalism based 
on the Nakajima-Zwanzig equation. Later time-local theories based on the same trick have 
been developed [10, 11]. For molecular wires these two versions of quantum master equations 
have been compared in Ref. [12]. 

Grossmann et al. [13, 14] showed in their pioneering work that periodically driven quantum 
systems can exhibit the unusual effect of quenching the tunneling dynamics at specific values 
of field amplitude. This effect was termed coherent destruction of tunneling (CDT). The 
validity range of this localization of particles by external periodic fields has later been studied 
in more detail for situations including a bias [15]. For cw-laser fields CDT has already been 
described for model molecular wires using the Floquet theory [16]. The purpose of this work 
is to study the CDT for short Gaussian laser pulses which have a high carrier frequency. 
Earlier studies [7] suggest that the CDT phenomenon survives also for shorter pulse length 
not just infinite length as in the original CDT studies. In those studies an average current was 
investigated. For cw-laser fields the CDT appears to be perfect if one measures at intervals 
separated by the period of the driving field but there are still oscillations within one period. 
Another way of looking at the effect is by averaging over one period of the periodic driving 
in the cw case. As will be discussed later, for shorter laser pulses one needs to average over 
several of the high-frequency oscillations of the carrier frequency. 

Model. - As in many quantum dissipation theories, the physical system is separated into 
the relevant system Hs(t), mimicking the wire, and reservoirs Hr modeling the leads 

H(t) = H s (t)+H R + H S R (1) 

with wire-lead coupling Hsr. The wire consists of N sites coupled to each other by a hopping 
element A (see Fig. [1]). Denoting the creation (annihilation) operator by c\ (c„) the tight- 
binding description of the electrons in the molecular wire reads 

H s( f ) =X!^™ c ™ c "~ A X!( c « c ™+ 1+c ^i c ") + ^5Z c ™ c ™ c ™+i c «+i~A l£; W • ( 2 ) 

n n n 
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The first term describes the on-site energies, the second term the nearest-neighbor hopping 
and the third term electron correlation within the wire in a version for spinless particles. The 
fourth term gives the coupling between the wire and the laser field E(t) whereas for the dipolc 
operator we assume [4] 

v-^ 2n - N - 1 + 
V = e Zs Xn = e 2^ 2 n ° n ' ( ' 

n n 

The environment of the wire consists of two electronic leads that are modeled by two 
independent reservoirs of uncorrelatcd electrons in thermal equilibrium. For each lead, the 
Hamiltonian Hr is given by Hr = J2 q ^q^ q c q with c q and c q creating and annihilating an 
electron in the corresponding reservoir mode \q) with energy u q . h is set to unity throughout 
the paper. Due to the assumed thermal equilibrium of the electronic leads, the occupation 
expectation values of the reservoir modes are determined by {c q c q i) = np(ui q — Ep)5 qq > where 
rip is the Fermi function and Ep the Fermi energy. In further derivations we will only refer 
to the left lead but the formalism has to be applied to the right lead coupled to the last site 
./V of the wire as well. The coupling of the left electronic lead with the first site of the wire is 
given by 

2 

H SR = K ^ = E(vk + n*4 c i) ( 4 ) 

x— 1 q 

with <!>! = J2 q V q c q , $ 2 = J2 q VqCq-, Ki = c i 7 ^2 = c i , an d a wire-lead coupling strength V q 
for each reservoir mode. 

As one is normally not interested in the dynamics within the leads but only within the 
wire, a quantum master equation (QME) based on a second-order perturbation theory in the 
wire- lead coupling has been developed for the reduced density matrix of the wire ps{t) [7] 

1,2 

p s (t) = -iC s (t)ps(t) - Y^ K ^ A ^'(*) - A xx /(t)] (5) 

xx' 

with Cs(t) being the Liouville operator applying Hs and the auxiliary operators 

A xx ,(t)= fdt'C^it-t'Wsit^K^psit'), (6) 

J to 

A**'(*)= [ dt'C* x , x (t-t')U s (t,t')p s (t')K x , . (7) 



J to 

Here we employed the definitions Us(t,t') = T + cxp ^— i f*, dTCs(r)j and the reservoir cor- 
relation functions C xx i{t) = trx{e tHRt § x e~ lHRt <& x < pn} with the reservoir density matrix pR. 
Using this form of the QME has the advantages over the standard Redfield approach [17] that 
the memory terms are included and that the wire-laser coupling is still treated exactly within 
the dipole approximation. 

For the QME all the external properties of the fermionic reservoir are described by a 
single quantity, namely the spectral density Jr(oj) = ^ q it\V q \ 2 8(w — u) q ). The sum becomes 
a smooth function for a dense spectrum of the reservoir modes. Now we use the trick of a 
numerical decomposition of the spectral density [7, 9] 
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with real fitting parameters pk, and Tk- This decomposition is not restricted to a any 
shape of the spectral density and can therefore be applied to approximate complicated band 
structures. This enables one to avoid the assumption of the wide-band limit and to take 
influences of the band structure on the dissipative electron transfer between the wire and the 
leads fully into account. 

With the complex roots of the Fermi function and of function ([8]) , the application of the 
theorem of residues results in 

m ~ . m m+m 

G «w = E snr + ^)e- in *') - f E ■feM)*-** = E ( Q ) 



771, 771 777 - \- 777 ' 

G »(*> = E ilV- ^ Ei?)ein * *) - 1 E = E ( 10 ) 

fc=l fc h P k k=l 

with the abbreviations fit" = ^fe + i^k and fi^ = O/. — iTk and the Matsubara frequencies 
Vk = i 27r ^ +7r + .Bp. Rigorously, the sum over the Matsubara frequencies would be infinite 
but it can be truncated at a finite m' depending on the temperature of the system T and the 
spectral width of Jr(w). The pure exponential dependence of the correlation function on time 
allows one to derive a set of differential equations for the auxiliary density operators 

^-A k xx ,{t) = 4 x ,K x , Ps (t) - i[H s (t),A k xx ,(t)] +1 k x ,A xx ,(t), (ii) 



g-fiAt) = (4,^ Ps(t)K x , - ilH s (t),A k x At)} + {riS ^>(t) (I 2 ) 

with A xx ,(t) = YTk=r' AL'W and A xx ,(t) = !X + r' A xx ,(t). 

As detailed in Ref. [7] it is now easy to give an expression for the current. Using the 
electron number operator of the left lead with the summation performed over the reservoir 
degrees of freedom AT; = J2 q c l c q yields 

h(t) = e—tr{Nip s (t)} = -ietr{[Ni,H(t)]p s {t)} - 2eRe(tr s {c\A 12 (t) - c\A 12 {t)^ . 

(13) 

Numerical results. - In the following we set the value 0.1 eV for the tight-binding hopping 
parameter A. This parameter also constitutes the energy scale for the results shown. The 
system temperature is set to T — 0.25A/fcs sa 290 K. Though with the method described 
above one can model complex energy dependences in the wire-lead coupling, here we use for 
simplicity only one Lorentzian in the spectral density ([8]). The maximum strength is chosen to 
be 0.1 A. This is achieved by using pi = 0.644 eV, Ti = 5.44 eV, and fii = 0.544 eV. With the 
chosen energy settings, a time unit in the system corresponds to 0.66 fs. The resulting current 
unit can be extracted from Eq. fT3")) and is equal to a macroscopic value of l[e] = 2.43* 10 -4 A. 
The carrier frequency is set to 1 eV. 

In Fig. [2] the phenomenon of CDT is shown for a monochromatic laser field of amplitude 
A = 2.405 eV. This amplitude corresponds to a full CDT [7,16]. The wire consists of two 
sites with on-site energies Ei = E 2 =0.54 eV which are initially empty. The Fermi energy of 
the left lead is 1.0 eV and that of the right one is 0.0 eV. The top panel of Fig. [5] shows the 
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Fig. 3 



Fig. 2 - The scenario of CDT for a monochromatic laser field. The top panel shows the laser field, 
the middle one the current through the wire, and the lowest panel the averaged current. 



Fig. 3 - Same as in Fig. [2] but for a short laser pulse with length a = 50 fs. 



laser field while the middle panel displays the time-dependent current I(t) through the wire. 
To be able to see the effect more easily we also show an averaged time-dependent current / 
determined by averaging I(t) over three periods of the highly-oscillating carrier field. Since the 
laser is instantaneously turned on at t = 0, one sees some transient behavior during the first 
150 fs before the CDT fully sets in. From Fig. [3] one can deduce that the transient behavior 
is mainly due to the slow onset of the CDT and not due to filling the wire sites. In this figure 
no laser pulse is applied to the wire for the first 150 fs and one observes a transient behavior 
during the first 50 fs due to filling the wire sites. Then the laser pulse with a Gaussian shape 
E(t) = A exp (— (£ — T) 2 / (2cr 2 )) sin (uidt) with a = 50 fs and a maximum amplitude as for the 
monochromatic laser before sets in. For the short laser pulse a complete suppression of the 
average current is not reached at any time but a dip in the current is clearly visible. The 
shape of this indentation does not show a Gaussian form nor is it symmetric. The current 
decays smoothly towards a minimum value while at the end of the laser pulse the current goes 
back to its equilibrium value showing some oscillations. 

To investigate the dependence of the average current on the pulse length, the average 
current is displayed in Fig. [3] for five different pulse lengths. It is clearly visible that the 
minimum average current approaches zero for longer pulse lengths. For an infinite pulse 
length the CDT condition would be perfectly fulfilled and the average current would vanish 
as shown in Fig. [2] It is interesting that already for a = 5 fs one can see a clear dip in the 
average current although only very few oscillations of the carrier field are involved in this 
scenario which reach an amplitude close to the maximum for that pulse length. 

All the calculations above were performed for U/A — 2. To investigate the effect of electron 
correlation on the CDT we also did the calculations using a Gaussian laser pulse with a = 50 
fs for different U. The results in Fig. [5] show already an interesting transient behavior of the 
current while the laser pulse is still off. Without correlation the equilibrium current has the 
largest value. Turning on correlation the equilibrium value of the average current decreases 
almost to the same value for electron attraction and repulsion oi U/A = ±2 . For the larger 
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Fig. 4 Fig. 5 

Fig. 4 - The average current through a molecular wire for different pulse lengths. Other parameters 
are the same as in Fig. 3. 

Fig. 5 - The effect of electron correlation for a pulse with length a = 50 fs. Displayed are cases with 
no correlation, electron attraction and electron repulsion. 



value of U/A = ±5 the equilibrium value of the current is reduced much more for electron 
attraction than for repulsion. Interestingly enough, the minimum value of the current is nearly 
the same in all cases regardless of the strength and sign of correlation. So the degree of CDT 
depends on the length of the laser pulse but not on electron correlation. 

Finally we want to show that effects of laser pulses are also visible if the conditions for 
CDT are not fulfilled. Then in general a much weaker suppression of the current is visible. 




Fig. 6 - After a short transient behavior no net current would flow through the wire for Ef,i = 0.82 
eV, Ef,t = —0.82 eV, E\ = —0.82 eV and E2 = 1.36 eV. The laser pulse which is the same as in 
Fig. [3] shifts the levels of the wire in such a manner that a temporary current is possible. 
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But if one uses a configuration in which, due to energetic reasons, a current is not possible 
but an external bias is applied, then there are cases in which the laser can initiate a current 
through the wire. This phenomenon can be seen in Fig. [6] The parameters are set so that 
the high energy level of the second wire site blocks the current. This creation of a current 
by a laser pulse does not seem to follow a pattern as the CDT. Its physical basis is quite 
different: The laser pulse brings the levels of the wire to a position which energetically allows 
for population transfer between the sites and between leads and outermost sites. So this effect 
can simply be explained by energetic reasoning and no interference effects are needed for an 
interpretation. 

Summary and outlook. - The main point of this work is to show that the phenomenon 
of CDT exists in molecular wires also for short laser pulses and not only for monochromatic 
laser fields. Nevertheless the amplitude condition as for cw-laser fields plays an important role. 
If these conditions are not fulfilled only a very weak current suppression is achieved by the 
laser pulse. If the CDT condition is fulfilled the current can be suppressed almost completely 
already for short laser pulses. For longer pulses the current suppression gets more effective. 
Electron correlation was treated in a version for spinless fermions. Large correlation can have 
a rather large effect on the equilibrium current, larger for electron attraction than repulsion. 
Finally we showed that for parameter regions in which no equilibrium current occurs, a laser 
pulse can lead to a weak current caused by temporary shifting the energy levels of the wire 
into a position in which population transfer is possible. 

In conclusion we showed interesting possibilities how ultra-short laser pulses can influence 
the current through molecular wires. Though the present calculations were performed for a 
simple model system, the basic physics should survive in more complex systems. Here we 
have restricted ourselves to a wire of only two sites although similar results can be obtained 
for longer wires. Using the effect of CDT with short laser pulses allows for the construction of 
fast opto-electronic switches. This opens, of course, a huge variety of applications for similar 
scenarios and possible technical devices even though a number of experimental hurdles have 
to be passed before one can think of direct technical applications. 
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